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Abstract: 

As is well known, the exact evolution equation of the curvature perturbation 
plays a very important role in investigation of the inflation power spectrum of the 
flat universe. However, its corresponding exact extension for the non-flat universes 
has not yet been given out clearly. The interest in the non-flat, specially closed, 
universes is being aroused recently. The need of this extension is pressing. We start 
with most elementary physical consideration and obtain finally this exact evolution 
equation of the curvature perturbation for the non-flat universes, as well as the 
evolutionary controlling parameter and the exact expression of the variable mass 
in this equation. We approximately do a primitive and immature analysis on the 
power spectrum of non-flat universes. This analysis shows that this exact evolution 
equation of the curvature perturbation for the non-flat universes is very complicated, 
and we need to do a lot of numerical and analytic work for this new equation in 
future in order to judge whether the universe is flat or closed by comparison between 
theories and observations. 
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1 Introduction 

In recent years the cosmology acquires a dazzling development. Unexpectedly, an 
elusory ingredient in the eyes of physicists, i.e., a tiny non-zero vacuum energy or 
so-called "dark energy" , turns out to be a main component of our present universe. 
Some evidence for the once recondite inflation has likely emerged now. A nice 
observed curve about the power spectrum of the Cosmic Microwave Background 
Radiation (CMBR) is actually fine consistent with the theoretical prediction [1]. The 
best model may be the slow-roll inflation of the flat universe. All seem faultless. A 
new era of the precision cosmology has been coming. 
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However, there is still some little bother. It seems the central value of the total 
density of the universe is slightly larger, £l T =1.02 ± 0.02. The large angle power 
spectrum looks like rather low [2] [3]. Whether or not these clues come into real 
questions? Maybe these are only observed errors, the future experiments will tell us 
that the universe is indeed flat, and the low power of the larger angle scales comes 
only from imprecise data. However, no matter how precise our future data are, we 
shall probably never be able to rule out the non-flat universe models. A couple of 
years ago, when one didn't find out enough matter density in the universe and didn't 
believe the cosmological constant will play important role in it, the open universe 
models were once voguish. Recently, the possibility of the closed universe model 
is being reconsidered [4]. Indeed, the idea of the closed universe model has a long 
history [5] [6] [7], some reasons seem attractive [8]. Anyhow, no matter the universe 
is flat, open or closed, it is important for us that we must know how to calculate 
exactly its primary power spectrum, in order for us to do a rigour comparison 
between observations and theories. 

As is well known, the inflation spectrum for the flat universe with single inflaton 
field has been understood very well, which is expressed by an evolution equation of 
the perturbation 

Xk + ( k2 ~ 4/ z o)Xk = 0, z = a<j)/H, (1) 
as well as the curvature spectrum 

k 3 lv, I 2 

and forms a standard paradigm[9], where the variable mass m 2 = —z'q/zq has an 
exact expression even for a no slow-roll case 

4/zq = a 2 H 2 {2 + 2ei - 3e 2 + 2e\ - 4 ei e 2 + e 3 ). (3) 

All notations can be found later. The fundamentality of all these results is never 
overestimated in the analysis of the primary power spectrum. 

However, for the non-flat case, especially for closed universe, the corresponding 
exact equations have not been written down clearly up to now, although some results 
are close to this goal nearly. Only when we know how to extend these equations to 
the non-flat case, and make theoretical predictions based on these new equations, 
we are able to judge whether the universe was flat, open or closed by observational 
confirmations. 

The goal of this letter is to provide the exact evolution equation of the scalar 
perturbation for non-flat universes from some elementary physical consideration. 
We do further some primitive, certainly immature, approximative analysis based on 
our exact equations. This letter is organized as follows. In section 2, we shall start 
with the perturbative formalism of the Einstein equation to obtain the curvature 
perturbative evolution equation for non-flat universes. In the section 3, we try to 
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do an analogic slow-roll analysis for these equations, look into whether the slow- 
roll approximation is suitable. In the section 4, we want to do a brief conclusion. 
Since our result may lay a foundation of further analysis for non-flat universes in 
future, the correctness of our result is very key. So it is necessary to give out the 
middle derivation as detailed as possible in order for interested someone to check it 
conveniently. In this paper, we only concentrate on derivation of the main formal 
results. A detailed analysis of the perturbation spectrum in phenomenology will be 
the subject of a forthcoming work. 



2 Perturbative evolution equation for non-flat uni- 
verses 

Our starting point is the standard effective action of the simple coupling system of 
Einstein gravity and single inflaton scalar field with an arbitrary inflation potential 

s = J Val-j^R + \<T%M,4> - v(<p)]d 4 x. (4) 

We shall take 8nG = M~ t 2 = 1. Now, supposing that inflaton field owns perturba- 
tion, we have 

<f>(t,x) = <f>o(t) +6<f>(t,x). (5) 

The perturbative metric in the longitudinal gauge, in terms of the Newtonian grav- 
itational potential $, is 

ds 2 = [1 + 2$(t, x )]dt 2 - [1 - 2<D(t, x)] °' (t) 2 (dr 2 + r 2 d6 2 + r 2 sin W)- (6) 

yL + K. C T /4J 

The background motion equations are 

# 2 + ^T = ^o + n 0o + 3F0o + ^ = O, H=-, (7) 

hereafter a "dot" denotes derivative with respect to the physical time, / = df/dt, 
and so on. The parameter K c takes 1, or —1 for the closed, flat or open universe 
respectively. 

From the standard literatures, we can get the Einstein tensor and the energy- 
momentum tensor in the perturbative form 

SG° t = + 5G° = - 3#<i> + 3(^ - H 2 )$, (8) 

(X (X 

ST? = i, ^ = - <£>] - (9) 

00 00 00 

Thus the linearized Einstein equation 5G\ = 4irG5T- for the component of Oi is, 



Another Einstein equation for the component of 00, having used (10), is 

A. = (1 + 2(™) )$ 

Both equations have been given in Ref.[10]. 

The motion equation of $ is simple. In order to express this equation in a more 
convenient way for later use we introduce a new variable S = a§. Using (10) and 
(11) we have an equation with regard to <3> in fact 

S = + a- 2 (A + 3K C )]S + j t [Ma<Pl)} • S. (12) 

However, what we need really is the motion equation about the intrinsic curvature 
perturbation R of comoving hypersurfaces, which can expressed as, seeing in Ref. [11] , 

R = $ + H^. (13) 

00 

In order to gain the exact equation about R in a clear way, we need to define the 
rolling parameters, as we emphasize, which are not only valid for slow-roll case, but 
also for no slow-roll case 

_ -2rr-2 _ 00 _ 00 _ 00^ _ 00^ n ,N 

„ = « ' « = 5h5. £2 = -^' c ' = Jpf (14 > 

The background motion equations (7) can be rewritten as 

H 2 = j(l + K cV - IcO" 1 , 0o = -^(1 - ^a)" 1 , (15) 

which formalism is particularly suitable for parameter expansion if the rolling pa- 
rameters are small, in order to be used later. During the progress of getting the 
equation about R, we shall meet many derivatives of these various rolling parame- 
ters with respect to time, which should be listed here, seeing Ref. [12], 

H — H 2 (K c r) — ei), r 1 = -2Hr ] (l+K c r ] -6 1 ), e ± = 2He 1 {e 1 -e 2 -K c rj), (16) 

e 2 = -H(e 3 -e 1 e 2 -e 2 2 + K c e 2 r ] ), e 3 = H[e 4 + e 3 (2ei + e 2 - 2K c rj)\. (17) 
Then the Eq.(12) can be rewritten as 

S = H 2 ^ + r](-k 2 + AK C )]S + H(l - 2e 2 )S, (18) 

where we have use the eigenfunction expansion for some fluctuation F(t,x) 

F(t,x) = J2 F ki m (t)Qki(r)Y lm (0,v), (19) 

klm 

(A + k 2 - K c )Q kl = 0. (20) 
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All definitions here can be found in Ref.[13]. Strictly speaking, k = \/k 2 — K c is just 
the dimensionless comoving wavenumber of a fluctuation. Hereafter our notations S 
and R mean their expanded amplitudes Skit) and Rk{t), since the isotropy implies 
independent on the indices I and m. 

Using Eq.(lO), the equation of curvature perturbation, i.e. Eq.(13), is recast as 

R = a - 1 S +(aHe 1 )- 1 S. (21) 

Doing derivative of it and using (18), we get 

R = ?^(ae 1 ) _1 [-(A; 2 - 4K C )HS + K C S}. (22) 

Solve R and R from above two equations, we have 

S = ae x {k 2 - 4K C + K^y^KM - aHR], (23) 

S = aHe x (k 2 - AK C + K^Y^k 2 - 4K C )R + aHe^}. (24) 
We can continue to apply derivative to (22), again using (18), yield 

R = (ah^ilK^ + (k 2 - 4K C )(3 + e x - 2e 2 - 2K c r])]S - H' 1 ^ 2 - 2K C )S}. (25) 

In the following we want to use the conformal time r = / dt/a, and hereafter a 
"prime" denotes derivative with respect to the conformal time, /' = df/dr, and so 
on. We have a transformation concerned R 

R = a^R', R = a- 2 R" - a^HR'. (26) 

Substituting Eqs.(23) and (24) into the right hand side of Eq.(25), and changing to 
the conformal time, we establish 

R" = -k 2 R + C R + dR', (27) 

where the coefficients are obtained finally after some calculation 

C = 5K C + 2K c {e l - e 2 - K c rj)[l + K c e x j{k 2 - AK C )\~\ (28) 

d = -2aH{K c e x + (1 + e x - e 2 - K cV )[l + K c e x j{k 2 - 4K.)}- 1 }. (29) 

The equation (27) has the first order term R', we hope to cancel it by introducing 
a new variable Xk, 

Xk = zR. (30) 

Thus equation (27) becomes 

Xl + k 2 Xk ~ C oXk - z' x z" Xk = 0, (31) 

if we take 

Ci = -2z^z'. (32) 
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Eq.(31) is a standard equation without the damping term and with an equivalent 
varying mass term m 2 = —C — z"/z. Solving z from (32) and (29), we get 

z = -a^e^l + K^/ik 2 - AK C )]-K (33) 

This is important parameter which controls evolution of the scalar perturbation. Its 
first and second derivatives with respect to the conformal time are given by 

z'/(aHz) = K c t x + (1 + ei - e 2 - K cV )[l + K c t x j{k 2 - ±K C )\~\ (34) 

z"/z = a 2 H 2 (k 2 - 4K C + K c e l )- 2 [k\2 + 2e 1 - 3e 2 + 2t\ - 4 ei e 2 + e 3 ) 
+K c k 2 (-16 - 12ei + 24e 2 - 15e 2 + 29eie 2 - 8e 3 - e? + 2e 2 e 2 - 3eie 2 + ei e 3 ) 
+if 2 (32 + I6ei - 48e 2 + 30e? - 52eie 2 + 16e 3 + 3e? - 8e 2 e 2 + 12e ie ^ - 4eie 3 ) 
+K c r](-6Se 1 - 4A; 4 ei + 32e 2 + 2A; 4 e 2 - 7ej + 16eie 2 ) + fsT c A; 2 ^ 2 (-16 - 
+K 2 k 2 r ] (33e 1 - 16e 2 + 2ej - 4 ei e 2 ) + X C V(32 + 2A; 4 + 4ei)]. 

(35) 

The equation (27) or (31) becomes finally 

X'l + [k 2 - 3K C - 2K c z'/(aHz) - z"/z] Xk = 0, (36) 

with an equivalent varying mass [— 3K C — 2K c z' /(aHz) — z" /z] 1 / 2 . Its spectrum, 
referring [14] and [15], is 

pR _ k(k 2 -K c ) \ Xk \ 2 

If the universe is closed, the parameter k, almost to be dimensionless comoving 
wavenumber, takes integer value 3, 4, • • •, to infinity. It is notable that the results 
gotten until now is exact, especially for (34) and (35), i.e., there is no any approxi- 
mation to be taken. Of course a precondition of such statement is to recognize the 
validity of the first order approximation for the perturbative Einstein equation. If 
K c = 0, all return to standard result for flat case, i.e., (1), (2) and (3), specially 
z — > Zq. Eqs.(33)-(36) are our main results, which are the foundation for further 
analysis. In fact, these results were able to be obtained almost by Refs.[15] and 
[16], if they would like to finish their last steps. However, we adopt a rather dif- 
ferent method to obtain the exact evolution equation of the comoving curvature 
perturbation for non-flat case, this increases the reliability of the results. The exact 
evolution equation of the tensor perturbation for non-flat universe will be found in 
a forthcoming work. 



3 A try to treat new evolution equation 

The evolution equation obtained by us for non-flat case is obviously more com- 
plicated than flat one. New parameter rj appears in an unmanageable way. How to 
treat it brings us a new challenge. What is a boundary condition on this equation? 
How to define the vacuum for this quantized system? Noted that rj = a~ 2 H~ 2 ~ r 2 , 
the coefficient of Xk term, i.e., variable mass square, not only contains a single term 
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of t~ 2 like the flat case, but also contains other terms which is very complicate func- 
tion of the conformal time r and even various e varied slowly. It is difficult to apply 
the slow-roll method in this case. We must think out some better methods to treat 
it. However, before finding out better outlet, we still hope to do a comparison if we 
father upon so-called slow-roll approximation on it. In such way we can appraise 
how the slow-roll approximation is partly invalid. 

Let us introduce the potential parameters which are very suitable for slow-roll 
approximation, seeing Ref.[12], 

Ul = (VJV) 2 , Vl = V H /V, v 2 = V^V^/V 2 , (38) 

to express the rolling parameters mentioned in (14), 

ei ~ \u x - \u\ + \u x v x + y K cUiV + |u? - |uf«i + + lu ± v 2 , , 

e 2 ~ \u x -vi- \K c r) - \u\ + |uiui - \v\ - \v 2 + 3K c mr] - fK^rj - §K 2 7] 2 
+§«? - 4u?ui + fu lV 2 - \v\ + y|miv 2 - \vivi - \v z - ^K c u\r) 
+ 7 -§K cUl v lV - K c v 2 r) - lK c v 2V + ^§K 2 Ul r] 2 - §K 2 v lV 2 + ^Krf, 

(40) 

e 3 ~ 2K c r] + u\ - \u x vi +v\ + v 2 - ^K^r] + y K cViV 
-fuf + fujvi - f u x v\ + \v\ - \u x v 2 + 2v x v 2 + |v 3 + fK c ujri - ^fK^^r] 
+ §K c v 2 lV + lK c v 2V + §K 2 V 2 - ^fK 2 c u x if + fK 2 v lV 2 + f K c r,\ 

(41) 

the terms with K c are our new results. Then we have, up to the second order, 



z'/(aHz) ~ 1 + ^ + -K cV - (-K c u lVl + -K^/ik 2 - AK C ) 

1 22 
+-(3u 2 - 9ui«i + 3v\ + 3v 2 - 13K cUl r] + 20K cVl rj + —K 2 r] 2 ), 
y o 

z"/z ~ a 2 H 2 [{2 + 7K c rj - \ Ul + 3^) - {\k c u iVi + 2i^ c 2 w 1?7 )/(A; 2 - AK C ) 
1 

+ — (33m? - 69miMi + 3Qv 2 + 36w 2 - U2K cUl r) + 150K cVl r] + 64K 2 r) 2 )}. 
18 



(42) 



(43) 



Moreover for the conformal time we have an approximation if the rolling parameters 
are small 

r ~ -a^H-^l + (ei - ^K cV ) + {3e\ - 2e x e 2 - ^K c e lV + ^K 2 rj 2 )}, (44) 

this result can be checked by dr/dt = a" 1 up to the appropriate orders by using 
Eqs.(16) and (17). If we take the slow-roll approximation, we have 

:(1 + o M i - T^ M i + - + —K cUl ri), (45) 



a# aiT 2 12 1 3 3 ' 9 
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then we get 



z"/z ~ r- 2 [(2 + -K cV + - Wl + 3 Vl ) - (-K c u lVl + 2K 2 u lV )/(k 2 - AK C ) 

3 9 3 2 2 19 4 ( 46 ) 

+ (o M ? + o M i^i + 2v l + 2 ^2 + 9K cUl r] + — A> 1? 7 - -rzK 2 c rf)]. 

The equation (36) is approximated as 

32 1 3 3 9 

Xk + - Y Kc + Kc ' °( Ul ' u< ) " ^t 2 + 2 Ml + 3ui + 2^ + 2 MlWl + 2u * 
3X^1^1 . . 9 

(47) 

+2v 2 - 2(^2-4^ ) J + Kc ' ( or ni S ner terms) }xfc = 0. 

Notice that we shall omit the 0(ui,Vi) terms in the term of r° order, and omit the 
higher terms such as r 2 , since they should be small in the limit case \r\ — > 0. Thus 
we obtained an equation very similar with flat case. Since we omitted the r 2 and 
higher terms, it is obvious that we are not able to do analysis on the limit case of 
\t\ — > oo. We then see that an Bessel asymptotic analysis loses its validity. 

If a Bessel analysis is forcedly adopted for the equation (47), then we have the 
spectrum 



where 



and 



32 

f(k) = k(k 2 - K c )(k 2 - -K c y 3 / 2 [1 + K^/ik 2 - AK c )l (49) 



3 1 5 7 12 1 

P = - + -u 1 + v 1 + —uj + -u lVl + -v\ + -v 2 - -K c u x v x j{k 2 - AK C ). (50) 

It is unnecessary for us to go further, not due to mathematical difficulty, but since 
the evolution is not slow-roll in the small k region, where is of greatest interest for us. 
In other words, we don't reassure our approximation, in spite of it is very efficient 
for flat case. However, these equations still give us some elicitation. We have seen 
that there will be some change in the forthcoming formulae of the spectral index 
and its running, which will be different from the flat case. Even if the universe is 
slight curved, the slow-rolling is still a good approximation for the power spectrum 
of the small angle part. 

We must note the relationship between the dimensionless comoving wavenumber 
k and the multipole I in CMBR spectrum if the universe is closed. If the horizon of 
our universe is about the present Hubble distance -f^o" 1 , the length corresponded to 
the multipole / is about ttHq 1 //, therefore this length should be equal to the physical 
half-wavelength na^jk corresponded to wavenumber k, thus we have k ~ clqHqI, 
where ao is the present cosmic scale factor with length dimensional. If the total 
density Vl T is about 1.01, then a H ~ 10. For the quadrupole fluctuation of / = 2 
in CMBR, which has the real cosmological signification, the wavenumber k is about 
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20. Therefore the enhanced effect of the factor f(k) for spectrum is not too large. 
The leading behavior of the spectrum is controlled by the parameter z ~ zq, 

k 3 x 2 - »T\ P? - ^V 1 - H±/4>1 (51) 

i.e., the factor H 4 /<pl in Eq.(48) governs the spectrum, just the same as the flat case. 
This key feature has been captured by Ref.[17], specially seeing their figure 15. They 
think this can explain the low power in small multipole / for CMBR spectrum (figure 
20 in [17]). 

If the universe is born from nothing to a finite size one through a tunnelling 
effect, as described by a model of Ref.[18], at this beginning the universe has an 
equilibrium of K c a~ 2 = V/3 from Eq.(7), we then have a asymptotic behavior near 
the time starting point t — > 

H ~ ~ t, H 4 /<p 2 ~ t 2 . (52) 

The quantity H 4 /<f) 2 will increase from zero up to its maximum at some time point, 
keeping this high value for a period with a very slowly decreasing, i.e., an approxi- 
mated de Sitter phase, then fall off to end its inflation. Before this time point the 
spectrum is suddenly low for small k, i.e., for a sky large angle or low multipole 
I, looks like a cut-off. It seems possible that the closed universe can explain the 
large angle lack of the CMBR spectrum and the little up-departure of the cosmic 
total density from the exact flat case. Anyhow we must remember that the result 
(47)- (50) obtained by us is only a try. In order to obtain an assured consequence, 
it is necessary to investigate numerically on the exact evolution equation (33)-(37), 
which will be our further serious tasks. 

4 Conclusion 

We obtained an exact evolution equation of the scalar perturbation for non-flat 
universe by using a direct method. The exact evolution equation of the scalar 
perturbation for flat universe has been playing an important role in the getting a 
primordial inflation spectrum. Its non-flat extension will play a key role in future 
to analysis whether the universe is closed or flat. The idea of the closed universe is 
attractive. After all, this idea is the simplest one among the many ideas to introduce 
a new scale. We should start with our investigation from as simple idea as possible. 

The exact evolution equation obtained by us is rather complicated. We have to 
do a try, in spite of which is immature, however, surely beneficial. It is impossible 
that the slow-roll approximation will lost wholly its efficiency even in non-flat case. 
Maybe some revelation can be drawn out from this rather simplified treatment. 
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